THE WEAK-TYPE (1,1) OF FOURIER INTEGRAL OPERATORS 

OF ORDER -(n-l)/2 



TERENCE TAO 



Abstract. Let T be a Fourier integral operator on R™ of order — (n — l)/2. 
In ji| it was shown (among other things) that T maps the Hardy space H 1 to 
L . In this note we show that T is also of weak- type (1,1). The main ideas are 
a decomposition of T into non-degenerate and degenerate components, and a 
factorization of the non-degenerate portion. 



1. Introduction 

This note is concerned with the mapping properties of Fourier integral operators. 
As our considerations will be local, we will work on R" for some n > 2, and we will 
define (as in Q ) a Fourier integral operator T of order m to be any linear operator 
of the form 

Tf(x) := f e™*^a(x,Of(0 « (1) 

where / is the Fourier transform 

/(0 := / e- 2 ™<f(x) dx, 

of /, a(x, £) is a standard symbol of order m (i.e., we have the bounds 
|^a(a: > OI<C a , /s (l + |e|r-l' J l 

for all multi-indices a, 0) with compact support in x, and $ is a real phase function 
which is homogeneous of degree 1 in £, is smooth in (x, £) for ^ on an open 
neighbourhood of the support of a, and obeys the non-degeneracy condition 

on an open neighbourhood of the support of a. 

A model case of Fourier integral operators arises in the translation-invariant 
setting, when &(x, £) = x ■ £ 4- (/>{£,) for some real smooth <j>, homogeneous of degree 
1, and when a(x, £) = a(^) is independent^ of x. Then T becomes a Fourier 
multiplier: 

Tf(0 - e 2 ^a(Of(0- (3) 

It is well-known (see e.g. M) that Fourier integral operators of order m < 
are bounded on L 2 . If m < — (n — l)/2, then one can also show 0] that these 
operators are bounded on L 1 and also on L°°. At the endpoint m = — (n — l)/2, 
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x Of course, this means that a is no longer compactly supported in x, but by applying cutoff 
functions as necessary we may make this operator consistent with the previous definition. 
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it was shown by Seeger, Sogge, and Stein 0] that Fourier integral operators map 
the Hardy space H 1 to L 1 . From this and Fefferman-Stein interpolation one can 
show |i) that a Fourier integral operator of order m is bounded on L p provided that 
1 < p < oo, || — 1| < —to/ (n— 1), except at the endpoint case when m = — (ra— 1)/2 
and p = l,oo. These results are sharpf]; see 01. 

This leaves open the question of what happens to operators T of order m = 
— (n — l)/2 on I 1 . One cannot expect T to be bounded on L 1 ; for instance in the 
model case (||), if = |£| and a is a generic symbol of order — (n — l)/2, then 
the convolution kernel if (a;) of T has magnitude comparable to l/(|x| — 1), which 
has a logarithmic divergence. However, we have 

Theorem 1.1. If T is a Fourier integral operator of order to = — (n — l)/2, then 
T is of weak-type (1, 1). 

The rest of the paper will be devoted to the proof of this theorem, but first we 
begin an informal discussion. 

We first consider the translation invariant case ([?]). To begin with we make 
the non-degeneracy assumption that V 2 </>(£) has the maximal rank of n — 1 for all 
£ 7^ 0. For sake of exposition we shall restrict ourselves to the case </>(£) = |£|, 
which is already typical. The principle of stationary phase then suggests that the 
convolution kernel K(x) of T should be concentrated near the image of V0, which 
in this case is the unit sphere. 

Let da be the standard surface measure on the unit sphere. As is well-known we 
have the Fourier transform expansion 

da(0=e 2 ^a + (0 + e- 2 ^a.(0, 

where a+, a_ are symbols of order — (n — l)/2 which behave asymptotically like 
a ±(0 ~ c±|£| _( - Tl ~ 1 * ) / 2 as |£| — > oo for some non-zero constants c±. If we ignore the 
e -27i-i|£| a ^ t erm (which can be suppressed by an appropriate use of cutoff functions 
in both space and frequency), we can thus factorize 

e 2 ^a(0=M0K0 
for some symbol b of order 0. This gives us a factorization 

T = SA 

where S is a pseudo-differential operator of order and A is the averaging operator 

Af := f * da. 

Since the A operator is clearly bounded on L , and since pseudo-differential op- 
erators of order are always of weak-type (1,1) by standard Calderon-Zygmund 
theory, we thus see that T is of weak-type (1,1) as desired. 

This simple example suggests that one should try in general to factorize the 
Fourier integral T into a pseudo-differential operator of order 0, and an operator 
bounded on L 1 . In the translation invariant case (||), this idea works well if the 
matrix V 2 $ is always non-degenerate (we will explain what uj is in the next section). 
However if V 2 $ degenerates, then the corresponding measure da develops cusps 
and it becomes difficult to control the asymptotics of do~(£). In principle this can 

2 More precisely, these results are sharp if V 2 </< has the maximal rank of n — 1, otherwise one 
can increase the range of p and m somewhat. Indeed, it is this gain in the degenerate case which 
allows us to prove Theorem We thank Michael Ruzhansky for pointing out this subtlety. 
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be avoided by damping the measure by a factor such as det(V 2 <F) 1/2 (cf. §, f§, 
I2J), however good asymptotics can only be obtained when |£| is very large (e.g. if 
|£| > (V**)" 1 / 6 for some small e). 

Thus it remains to handle the contribution when <j> is close to degenerate, and 
when |£| is not too large. However, it turns out that this portion of the Fourier inte- 
gral operator is in fact bounded on L . For instance, in the completely degenerate 
case when $(£) = xq ■ £ is linear in £, the convolution kernel K is essentially a frac- 
tional integral kernel l/\x — xq\, which is clearly integrable. More generally, when 
<& is close to degenerate, then the error terms in a Taylor expansion of <& become 
more favorable, and one can coarsen the standard "second dyadic decomposition" 
(see e.g. |], §) in order to improve the standard estimates on the kernel K 
(which in the non-degenerate case, just barely fail by a logarithm to be integrable, 
because the Fourier integral operator has the critical order — (n — l)/2). 

To summarize, our strategy in dealing with a general Fourier integral operator 
T will be to first decompose T into a "degenerate part" (roughly, this is where 
l£l det(V 2 <I>)~ 1 / e ), and a "non-degenerate part". The kernel of the degenerate 
part can be shown to be integrable, which is acceptable. The non-degenerate part 
will be factored (modulo errors which are essentially Fourier integral operators of 
order strictly less than — (n — l)/2) as the product of a pseudo-differential operator 
of order 0, and an operator whose kernel K(x,y) is essentially a (smoothed out) 
measure on the singular set {(x, V{<i>(x, £)) : (x,£) G supp(a)} corresponding to 
<&, damped by the natural factor of det„_i (V 2 ^) 1 / 2 . In particular, the kernel K 
is integrable. (The smoothing out will arise from the cutoff to the region |£| ^> 
det„_i(V 2 <£ > )~ 1 / £ , but this will not affect the integrability of the kernel). 

This project initially began in 1995 while the author was a graduate student at 
Princeton, and was suggested by the author's advisor, Elias Stein, who in particular 
suggested the use of factorization. The work was finally completed at UCLA and 
ANU. The author is deeply indebted to Elias Stein for all the support, encourage- 
ment, and guidance given over the years. The author also thanks the anonymous 
referee and Michael Ruzhansky for helpful comments. The author is currently a 
Clay Prize Fellow and is supported by the Packard Foundation. 

2. Notation and preliminary reductions 

The dimension n will always be fixed. We shall need a small number < e <C 1 
depending only on n; for instance s := n~ 10 will suffice. When a constant C appears 
in front of e, e.g. 2 Csk , it is understood that C is independent of e and depends only 
on n. (In particular, Ce can be made arbitrarily small by choosing e sufficiently 
small). 

We use T*R" := {(x, £) : x, £ £ R"} to denote the cotangent bundle of R™. We 
shall use the Euclidean metric and standard basis e±, . . . , e„ throughout, and so we 
will not bother very much to distinguish between a vector and covector, or between 
a matrix and a quadratic form, etc. 

We will use A < B or A = O(B) to denote the estimate \A\ < CB, where C 
is a constant depending only on n, e, the C K norms of $ on the cosphere bundle 
{(x,£) G T*R" : x,£ G R n ;|£| = 1} (where K = 0(1) depends only on the 
dimension) and the constants in the symbol bounds for a. 

We will write the frequency variable as £ = (£, £ n ), where £ G R™ -1 and £„ G 
R. We similarly decompose the spatial variable x = (x,x n ). We shall also use 
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"projective co-ordinates" (A,w), denned by 

A = £«; w = £/£„, 

to analyze the frequency variable £ (if £ is close to the e n axis). Even though A and 
are equal in value, we shall distinguish between the radial derivative d\ (which 
keeps u> fixed) and the vertical derivative d^ n (which keeps £ fixed) . One of the key 
features of the argument will be that all the cutoff functions are very smooth with 
respect to the radial derivative d\ even if they become rather rough in the angular 
directions V w . 

Fix a, $. We now make some basic reductions to simplify the form of T. 

We first observe that we may assume that the symbol a(x,£) vanishes when 
£| < 1. This is because the portion of the operator on the region |£| < 1 is a Fourier 
integral operator of arbitrarily large negative order, and is therefore bounded on L 1 
by (for instance) the results in Q] . The main purpose of this reduction is to ensure 
that we will not have to worry about any possible singularity at £ = 0. 

Next, we can assume after a finite partition of unity of the frequency angular 
variable £/|£| and a rotation that the symbol a(x,£) is supported on the cone bundle 

C:={(x,£)eT*R n :£„»l;|£|«A}. 

This allows us to use A as a proxy for |£|, and u> as a proxy for the angular variable 
This will be convenient technically because the hyperplane R™ -1 is flatter 
than the sphere S n . From the homogeneity of the phase function $ we have 

$(a:,0 = A$(a:,(w,l)). 

We will abuse notation and write $(x,uj) for $(x, (w, 1)). 

We will always be working on the support of a(x, £), so we shall implicitly assume 
that (x,£) € C throughout the rest of the paper. 

Define the canonical relation E C (T*R n ) x (T*R™) by 

S := {((V e $(z, 0, 0. & V*$(ir, 0) : (x, G supp(a)}. 

We will not use S directly, but the geometry of this relation will be implicit in our 
arguments. Observe from (^) that by restricting the support of a if necessary we 
can make E an embedded manifold, and in particular we can assume that the map 
(ai, £) i— > (x, V x <l'(x,£)) is a diffeomorphism on the support of a. From (0) and the 
Euler homogeneity relation 

V x *(a:,0=e-V c V B $(a:,0 
we observe the useful bound 

iv, $(£,£) i (4) 

Fix x £ R n and u> & R n_ . A key quantity in the analysis will be the curvature 

J{x,u) :=det(V£$(a;,w)). 

This curvature measures the extent to which the phase function fails to be linear; 
thus $ degenerates when J is small. Observe that J is smooth in x and u>. 

We shall frequently need smooth dyadic cutoffs of Littlewood-Paley type. We fix 
</>(£) = </>o(£) to be a non- negative radial bump function on R" which is supported 
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on the ball {|£| < 2} and equals 1 on the ball {|£| < 1}. For any k <E R, we define 
the functions (f>k, i]k by 

MO ■= 0(e/2 fc ) 

%(0 :=0fc(O~0fc-i(O- 



By abuse of notation we shall also define these functions on R" 1 and R in the 
obvious manner. 

3. Splitting into degenerate and non-degenerate components 

Let T be as in the previous section. We can decompose T into Littlewood-Paley 
components^] 

where the restriction A; 3> 1 is justified since a(x,£,) is supported on the region 
We can thus split T = Td cg + T non ^ eg where 

T dcg f(x) := ]T f e 2 ^^a(x,C)cf>- ek (J(x,oj))r, k (Of(0 d£ 
fe>i 

and 

T nondeg /(z) := L e* ri *t x '«a( a:> 0(l - ^(A*. C 

fc^i- 711 (5) 

Informally Td og is the portion of T where |J(x,o;)| < and conversely for 

?nondeg- Thus Td og captures the portions of T which are even just a little degenerate. 

We can now outline the remainder of the proof. In Section || we will show that 
Tdeg is bounded oni 1 : 

IIWIIi S Wfh- 

This will basically be a straightforward computation of the kernel of Xd og , modifying 
the standard second dyadic decomposition slightly. 

To deal with T non d eg requires more work. We will obtain a factorization 

Inondeg = SA + E 

where A, E are operators bounded on L 1 

||Aflli,l|£/lli<ll/lli 

and S is a standard pseudo-differential operator of order (so in particular is of 
weak-type (1, 1), see e.g. 0). From all the above estimates we immediately obtain 
that T is of weak- type (1,1). 

The operator A will have a kernel K(x, y) which is essentially a finite measure 
on the set 

(i.e., the spatial projection of the canonical relation £). The measure is weighted 
by the natural damping factor |</(y, £)l^ 2 , but with some additional cutoffs and 
blurring arising from the truncation to the region |J(x,o;)| 3> The operator 

A is thus the natural averaging operator corresponding to the canonical relation S. 
We construct A and prove its L 1 boundedness in Section 

3 All summations will be over the integers unless otherwise indicated. 
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The operator A will turn out to essentially be an Fourier integral operator of 
order — (n — l)/2 with the same phase function as T. From the symbol 

calculus of Fourier integral operators one can then solve for 5*, in such a way that the 
error E is essentially an Fourier integral operator of order 1/2 better than — (n— 1)/2 
(except that the symbol estimates have been slightly worsened by the presence of 
the cutoffs 4>- e k)- We perform this procedure and prove the L 1 boundedness of E 
in Section |6|. 

4. Boundedness of T deg 

We now prove the T^-boundedness of Td cg - By the triangle inequality it suffices 
to show that 

|| f e 2 ^^a(x,0<l>-sk(J(x,u;))Vk(Of(0 < 2^ 

for all h > 1. 

Fix k. By Minkowski's inequality it suffices to show that 

|| / ^^^(..^^(Jta;^))*^^ 1 ^ ^IUi <2- £fc (6) 

J XV 

for all y. 

Before we prove (^), let us first informally review the proof (from Q|) of the more 
standard estimate 

|| f e^^a(x,0vk(0e- 2nty< #11*1 < 1; (?) 
J R 

this estimate does not restrict the integration to the degenerate region, but on the 
other hand it does not obtain the crucial decay of 2~ ek . In other words, (0) asserts 
that the Littlewood-Paley pieces of T are uniformly bounded on L 1 . 

To prove (Q) we apply the "second dyadic decomposition" and partition the uj 
variable smoothly into about 2 (n ~ 1 ) fc / 2 disks D of radius 2~ k / 2 . It then suffices to 
show that 

|| [ e 2 ^ x ^a kD (x,0e~ 27Tlv< d^Li <2~ (n ~ 1)k/2 (8) 

for each of these disks B, where ak,D is the symbol a smoothly cut oflQ to the 
tubular region {(x,£) £ C : A ~ 2 k ;uj e D}. Note that the £ variable is now 
restricted to a tube of volume ~ 2(" +1 ) fe / 2 . 

We continue our informal discussion. Let ujd denote the center of the disk D. 
Observe from homogeneity and the Taylor expansion that 

*(x,O-I/-£ = A($0c,w)- !/•(",!)) 
= \{${x,w D ) - y (u D ,l) 

+ (V u *(x, cod) -y)-{uj~ ojd) + E(x, £) 
where the error term E(x,£) has a leading term of 

^VI$(x,oj d )(uj - uj d ,oj - oj d ) (9) 

4 Note that in applying this cutoff, a large portion of the angular regularity of the symbol a is 
destroyed, in that the angular derivatives are much larger when applied to £, rather than 
aj.. To put this another way, much of the angular regularity of is superfluous. This spare 
regularity in the angular directions will be crucial in all of our arguments. 
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(where we think of V 2< J> as a quadratic form). 

Because A = 0(2 k ) and to — ujo = 0(2~ fe / 2 ), the contribution of the error term 
E(x,£) to the phase Q(x,£) — y ■ £ is 0(1), which can then be absorbed into the 
amplitude function a^.D- From the principle of stationary phase we thus see that 
the integrand in (ph should be extremely small unless we have 



$(x,cj D )-y(uj D ,l) = 0(2- k ) 



and 



V u $(x,uj d ) - y = 0(2- fe / 2 ). 

Because of the non-degeneracy condition (||) , this restricts the variable s to a diskQ 
On this disk we can crudely estimate the integrand of flq) using absolute values by 



\~2 k ;ujeD 



the claim (||) (and hence (fy) then 



Since the volume of the disk is 0(2"(«+ 1 ) fc /2) 
follows. 

To improve (R) to (^J) we have to take advantage of the degeneracy condition 
J(x,lu) = 0(2~" : ). The point is that this degeneracy allows one to estimate the 
leading order error (^]) slightly better, especially if ui — ujjj lies in a direction where 
\7^&(x,ujd) degenerates. This allows us to widen the disks D used in the second 
dyadic decomposition to slightly larger ellipsoids while keeping the phase error (|^) 
of the order of 0(1). Intuitively, this should let us decompose into fewer tubular- 
regions, which will be the source of the 2~ ek gain. 

There is an apparent technical difficulty (when n > 3) in that the ellipsoid around 
u>d depends on the eigenvalues and eigenvectors of the symmetric n-lxn-1 
matrix V 2 $(x,wc), and so the eccentricity and orientation of the ellipsoid will 
vary with the center lod- This raises the fear that one would need a Kakeya-type 
covering lemma to partition phase space properly, but fortunately the above vari- 
ations are very smooth, so that nearby ellipsoids have almost identical orientation 
and eccentricity, and one can use a standard continuous partition of unity to obtain 
the desired estimate (^). 

We now turn to the details. We shall formalize the argument as the following 
Proposition: 

Proposition 4.1. Let k 3> 1, y € R™, and let afc(x, £) be any function supported 
on the region 

{(x,0 e supp(a) : A ~ 2 k , \J(x,u)\ < 2^ k } 
which obeys the symbol-type bounds 

\dgffla k (x,£)\ < c Pn 2-( n -^ k / 2 2-^ k 2 CEk W (10) 

for all multi-indices /3,7 (regularity in x will be unnecessary). Then we have 



R 



^i(H*,i)-yi) ak (x,0 



<2 



-ek 



(ii) 



5 This disk can be viewed as the projection of the above tube in phase space via the canonical 
relation S. 
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The estimate (^|) follows immediately from this proposition by setting 
a>k{x,Q := a(x,£,)(j>- £k (J(x,u}))rjk(Q. 

Observe that the cutoff <fi- e k causes the unavoidable 2 ek loss in the symbol es- 
timates (|hJ). These losses mean that we lose control of the regularity for scales 
of ui greater than 2~ Cek , but this will not be dangerous as all of our stationary 
phase computations will be on tubes for which the fluctuation in lj is at most 
0{2 Cek 2- k / 2 ). 



Proof (of Proposition 4.1) We first begin by defining a positive-definite analogue 
of the matrix V;J,<I>, which will be used to define the ellipsoids mentioned earlier. 
For all x,lu, let Q{x,u) denote the positive-definite u-lxn-1 matrix 

Q(x,u) := (2^ k ld+(Wmx,u J )) 2 ) 1 / 2 

where we now think of V^*!* as a real symmetric matrix, and Id is the n-lxn-1 
identity matrix. Observe that as quadratic forms we have the estimates 

2- fe |c| 2 + |JM(C,C)I < QM(CC) < ICI 2 , 

so in particular Q dominates J. We also have the smoothness property 

\d*Q\ < C P 2 C ^\ 
for all multi-indices (3. Finally, we observe that 

2- Cek < det Q(x, u) < | det J(x, u)\ + 2" efe . 
Let x € R™ and ujd € R™~ . We define the function ^ x ^ d (uj) by 
2(»-i)*/a0_ fc (Q( a . ) Ud )(uj -uj d ,lu- lu d )) 



det(Q(x,tj D )y/2 



where 4>-k was defined in Section B observe that this function is L^normalized 
and supported on some ellipsoid centered at u>d which contains the disk {uj — lud + 
0(2~ k / 2 )} but is contained inside the slightly larger disk {uj = uj D + 0{2 Cek 2~ k / 2 )}. 
We then define the averaged function i(j x (uj) by 



"4> x {uj) := I Vs,w D (w) duo; 

note that this is well defined for (x,^) in the support of afc and k sufficiently large, 
since we must have ujd = w + 0{2 Cek 2~ k / 2 ) by the previous discussion. 

The function tj) x is clearly positive. In fact we have the following estimates on 

Lemma 4.2. If(x,^) is in the support ofak, and k is sufficiently large, we have 

Ipx(uj) ~ 1 

and 

l^xHI <C p 2 Cek ^ 

for all multi-indices /3. 
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Proof We begin with the first estimate. We compute 



CgR" 1 :C=0(2C»fe2- fc / 2 ) 



2(n-l)fc/2 



det(Q(x, oj + 0)- 1/2 ^k(Q(x, oj + C)(C, 0) dC- 

'C=0(2°**>2-k/ 2 ) 

By Taylor expansion and the constraints on £ and Q{x,uj) we have 

det(Q(x, u + C)r 1/2 - det(Q(a; ) lo))~ 1/2 + 0(2 Cek 2- k ' 2 ) 

and 

^ k (Q(x, w + C)(C, 0) = 4>-k(Q(x, w)(C, 0) + 0(2 C£fc 2- fc / 2 ); 
It is clear (from the estimates on the non-degeneracy of Q) that the error terms 
0{2 Cek 2~ k / 2 ) will contribute at most 0(2 Cek 2- k / 2 ) to ip x (uj); this will be accept- 
able if e is sufficiently small. This leaves the main term 

2t™-W2 det (Q( a; , w ))-i/2 f (j)_ k (Q(x,uj)(C,0) dC 



But by a change of variables this is equal to J </>o(CC) d(, which is some positive 
absolute constant, and we are done. 

The second estimate is proven similarly to the first; when one applies an u> 
derivative, the main term J <^o(CjC) d£ disappears (since it is just an absolute 
constant), and the error term gets larger by about 2 Cek for each derivative. We 
omit the details. ■ 



We now return to the proof of (fLl|). We split the integrand as an average of 
integrals over ellipsoids: 



R 



R 



R" i>x(w) 



Since u>d effectively ranges over a compact set, it suffices by Minkowski's inequality 
to show that 



e 2m($(x,g)-yg) ak{x,Q4>x,u D (u) ^ 
R" i>x(u) 



<^ 2~ ek 



Li 



for all u> d (this is the analogue of (@)). 

Fix w£>; we may assume that Q(x,ui£>) = 0(2~ ek ) since the integrand vanishes 
otherwise. We split £ into polar co-ordinates A, u> and rewrite as 

ic^p J ^^^^ I 



We make the change of variables 

cj = cj D + 2- fc / 2 Q( a ;,^)- 1 / 2 C 
and expand out ?p x .uj D (u>) to rewrite the previous estimate as 



< 2 



-efc 



(12) 



<^ c^ — iri — X^k jlc^—ek 

Ll " ' (13) 
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where 

h (xn ,-(»-i) fc / 2 A"- 1 a fc (x,A(u;,l))0o(|C| 2 ) 

bk,UD,x(\ C) ■= 2 ^ 

and u> is of course now given by (|l^) . 

The amplitude bk,u D ,x is supported on the region {A - 2 fe ;C = 0(1)}. Because 
of the normalization factor 2~(™~ 1 ) fe / 2 and the estimates on a k , ip x we have the 
symbol bounds 

\d& s c b k , UDtX (\, C)| < Cp, s 2- k W 

for all /3, 8 (if e is sufhciently small), because any losses of 2 Cek in differentiation 
in the £ variable are always accompanied by gains of 2~ k / 2 thanks to the change of 
variables (|l2|). 

We now expand the phase in a Taylor series as before, to obtain 
w) - y ■ (uj, 1) =$(x, lo d ) - y ■ (wd, 1) 

+ 2- fe / 2 C • Q(x, w)- 1 / 2 (V w $(x, u D ) - y) 

(0 

where the error ek,x,y,u D (C) h as an expansion 

d (oj) = l2- k Vl^x,u JD )(Q(x,u JD )-^\,Q(x,u JD )-^\)+0(2 c ^2-^ 2 ). 

Since Q dominates V 2 $, it is easy to obtain the bounds 

\d s c e kiXt y !U)D (u;)\ <C s 2- k 

for all multi-indices S; once again, any losses of 2 Cek are compensated for by gains 
of 2~ fc / 2 . Thus the phase term exp(2ni\ek,x ,y ,u D ) can be harmlessly absorbed into 
the amplitude function bk,u D ,x- 

By the principle of non-stationary phase (see e.g. [Q), we can thus bound the 
integrand of ( |l3| ) pointwise by 

2 k {l+2 k mx^ D )-yiu D ^))\)- w ^\l+2 k l 2 \Q^ 

This function has maximum size 2 k , and is rapidly decaying outside of the eccentric 
disk 

{x : u D ) = y ■ (uj d , 1) + 0(2- fe ); V u $(x, = y + Q(x, cu D ) 1 / 2 0(2' k / 2 )}. 

Because det Q(x, w D ) < 2~ efe , this disk has volume 0(2- ek 2- (n+1 ^ k / 2 ). The claim 
(H) follows. ■ 

The proof of the L 1 boundedness of Td og is now complete. 

5. Construction and boundedness of A 

We now construct the averaging operator A. This operator needs to be bounded 
on L 1 , and also essentially be a Fourier integral operator with the same phase 
function as T. To motivate matters, let us first suppose we are in the non- 

degenerate case, so that \J(x, u>)\ ~ 1 throughout. Then the most natural choice 
for A is the operator Aq defined by 

A f(x) := J f(y&(x,w))<p(x,u) dio 

where V^$(x,w) is the function V^$(x,0 evaluated at £ = (uj, 1), and ip(x,uj) is a 
suitable bump function to be chosen later. One can motivate the choice of Aq from 
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an inspection of the canonical relation £, but we will instead use the principle of 
stationary phase in the informal argument which follows. 

It is clear from (0) that Aq is bounded on L 1 , since for each u> the map x i— » 
W^(x.lu) is a local diffeomorphism. Now let us write Aq as a Fourier integral 
operator. From the Fourier inversion formula we have 

In accordance with the principle of stationary phase, we now look at where the 
phase is stationary in lu: 

From homogeneity we have 

V f *(a:,w) = (V€*(x,(w,l)),^ n *(a!,(w,l))) 

= (y w ®(x, w),$(x, oS)-uj- V w $(a;, uj) 
and hence we have the identity 

V U K' • V c $(a;, w)) = V w (AV ■ V w $(x, w) + \'(&(x, u)-w V w $(x, w))) 

= AV-w)-V*$(a:,w), (14) 

where we of course write £' = (£', and A' = w' = Since we are in the 

non-degenerate re gim e \J\ ~ 1, we thus see that the only stationary point occurs 
when u) — lu'. By (p~4|), the Hessian at this stationary point is 

det(V5(£' • V^fow)))^ = det(-A'v50>(x, W ))|^ = {—\' n ) n ~ x J(x,w') 

(15) 

and the value of the phase at the stationary point is 

2iti£' • V 5 $(x, lu') = 2ni^(x, £') (16) 
so by the principle of stationary phase (see e.g. 0) we have the asymptotics 

= e 2^*(^«')^( XjU; ') e M^V)W4 A -(«-l)/2| w /)|-l/2 + ... 

for large where fi(x, lu') is an integer quantity which depends only on n and the 
signature of V^<I>(:r, u>) (and in particular is a constant in the non-degenerate case 
|J| ~ 1). Ignoring the error term (which will give a Fourier integral operator of 
lower order), we thus see that Aq is an Fourier integral operator of order — (n — 1)/2 
with phase £) and symbol given by 

^(x, W ')e^ /4 (A')- (n - 1)/2 |J(x, W ')r 1/2 . 

This choice of A will work in the non-degenerate case, but as J becomes degen- 
erate we see that the symbol of Ao can get large, so that Aq is no longer of order 
— (n — l)/2. The obvious remedy is to damp Aq when J is small, and use a new 
operator 



Aif(x) := / f(\7^( Xl u;)Mx,Lo)e-^^ m / 4 \J^^)\ 1/2 dcu 

2 ^'- v ^ x '^<p(x > u)e-^ x ' uj)ni/4 \J(x,u))\ 1/2 dcu) /(£') 



e 2 * 
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since this (heuristically at least) will be an Fourier integral operator with phase 
^(x,^) and amplitude 

ip{x 1 uj){X')- {n - 1 ^ 2 . 

(Compare this with |Q, ||). However, the integrand in the above definition of A\ 
has singularities on the zero set of J{x,u), and so we need to apply a cutoff away 
from that region, compatible with the cutoff (||) used to define T non( jeg- 

From the above informal discussion, it is now natural to construct A via the 
formula 

Af(x) := V [( / e 2 ^'- v ^^^V(^^)e" M ^" )W4 (l - <f>-ek{J(x,u>))) 



(compare this with (||) and the definition of A±). Observe that the cutoff (1 — 
<t>-ek(J{x,u>)) ensures that J] 1 / 2 and fi(x,uj) are smooth functions on the region 
of integration; indeed, the integer-valued function /i is constant on each connected 
component of this region. 

In the remainder of this section we show why A is still bounded on L 1 . In 
the next section we explain why we can factorize T non d C g = SA + E where S is a 
pseudo-differential operator of order and E is bounded on L 1 . 

We make the a priori assumption that / is smooth and compactly supported; this 
assumption can be removed by the usual limiting argument. We can also assume 
that / vanishes for |£| < 1, since on this region of frequency space, k is bounded 
and it is easy to show the i 1 -boundedness of A in this case. 

We write = </>*:(£') — ^k-iiC) an d apply summation by parts, to rewrite 

Af(x) as 



E /(/ e 2 ^'- V ^ x ^ip(x,u;)e-^ 
fc»i J J 



X,L))lti/4: 



(<l>-e(k+l)(J( x >u)) - <j>-ek{J{x,Uj))) 



(the vanishing properties of / ensures that there are no boundary terms). It will 
thus suffice to show the estimate 

|| JU e 2 ™«'- v «*^V(^)e~ M( ^ W4 

(4>-E(k+i)(J(x,u)) - <p- £k (J{x,uj))) 

<t> k {o\JM\ 1 ' 2 doj)na dehi z 2- £fc / 2 n/i| 1 , 

since the claim then follows by summing in k using the triangle inequality. 
Fix k. By the Fourier inversion formula we can rewrite the left-hand side as 

(4>-E(k+i)(J{x,u)) - cj)- £ k(J(x,uj))\J(x,uj)\ 1/2 duj\\ L i 
where Pkf is the Littlewood-Paley operator 
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Taking absolute values everywhere, we can estimate the previous crudely by 

/ / \P k f(V$Hx,w))\\(p(x,u)\2- ek / 2 dxdu. 
By (|^) one can bound this by 

f 2^ k / 2 \\P k f\\ 1 du ~ 2^ 2 \\P k f\\ 1 < 2- £fe / 2 ||/|| 1 , 

J\u\<l 

as desired. 

6. Construction of S, and boundedness of E 
We now need to find a suitable pseudo-differential operator 

Sf(x) : / e 2 ™<s(x,£)f(0 ,/£ 



where s satisfies the standard symbol estimates of order 

\dy^(x,o\<c a , a + \^r m (i7) 

for all multi-indices a, /3, such that we have a good factorization of the form 

Inondeg = SA + E. 

We begin, as usual, with a heuristic discussion. By repeating the stationary 
phase computations of the previous section we see that A should essentially be a 
"Fourier integral operator" with phase (fr(x,£,) and symbol 

^(x,u;)X-^/ 2 (l - 0_ efc (J(a:, w ))) % (O. 

We enclose "Fourier integral operator" in quotes, because the symbol does not 
quite obey standard symbol estimates (there is a slight 2 ek loss of regularity in the 
angular variable oj). From standard symbol calculus (or more stationary phase) 
we thus expect SA to also be a "Fourier integral operator" with phase <E>(a;,£) and 
symbol 

fc>l 

Meanwhile, T n ondcg is a "Fourier integral operator" with phase <i>(x,^) and symbol 

^2 a(x,0(l - <l>- e k(J{z,(j))rik(Q. 
fe»i 

Thus, if we choose <p to be bounded away from zero on the support of a, and define 
s implicitly by 

\(n-l)/2 I t\ 

s(x,V x $(x,Z)):=- *?l5L (18) 

(p{x,u) 

on the support of a, and define s to vanish otherwise, then we see (from (@) and 
the hypothesis that (x,£) — > (x,V x &(x,£,)) is a diffeomorphism on the support of 
a) that s does indeed obey the symbol estimates (17). 

It remains to show that the error operator E := T non deg — SA is bounded on L . 
It suffices to show that 

\\TnondcgSz - SAS z \\l < 1 

for all z e R" , where S z is the Dirac delta at z. 
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Fix z. We can expand SAS z (x) as 

^ SIS / e2 " [( ^ y) ' C+C '' (V5 * (y '" ) " Z)ls (^,C)^(2/^)e- Alfe '" )W4 
fc»i ■> J J J 

(1 - 0_ £fc (J(y, W )))77 fc (O|./(y,w)r /2 dudCdt'dy. 
On the other hand, from ( p"8| ) and (^) we can write T non degd z (x) as 

It thus suffices by the triangle inequality to show that 



(19) 

for all k ^> 1, where 

(l-^CJfe,^)))!^,^)! 1 / 2 d^cdy 

and the phase ^ x ,k.z is given by 

*, lM (u;, C, y) :=(z - 2/) ■ C + • V ? $(y, w) - <f(x, £') 

=(« - y) ■ C + r • (V 4 *(y, w) - V s $(x, £'))• 

The function W° is basically a symbol of order — (n — l)/2, modulo errors of 
0(2 Cek ). What is not as obvious is that the quantity W — W° is a symbol of slightly 
smaller order. More precisely, we have 

Proposition 6.1. We have the estimates 

\d^(W x . kiZ - W£ ki J(Ol < C /5 2 c ^ fc 2- fe / 2 2-(™- 1 ) fe / 2 2- fc l /3 l (20) 

when |£' ~ 2 k . 

In other words, apart from some errors of 2 Cek , W — W° is a symbol of order 
— (n — l)/2 — 1/2. Assuming this proposition, we can repeat the derivation of (0) 
(taking some losses of 2 Cek ) and thus bound the left-hand side of ( |T9| ) by 2 Cek 2~ k / 2 , 
which is acceptable if e is sufficiently small. Thus to complete the proof of Theorem 
1.1 we just need to verify the proposition. 



Proof (of Proposition 6.1) Fix x, k, z, We begin by using the principle of non- 
stationary phase in the u>, y, and C variables in turn to truncate the integrals in W 
substantially. 

We first observe from (|TJ) that 

Vu*x,k,z(u, C, y) = - u) • V 2 $(x, uo). 
Since the integrand is supported on the region where \J(y, > 2~ efc , we thus have 

|V w * x , M (w,C,2/)| > 2V -u\. 
Thus if we insert the cutoff 1 — 4'-(i/2-e)k{ UJ ' ~ w ) m the integrand of W, we 
have \Vu)^x,k,z{u, C>y)| <L 2 ek 2 k / 2 . Repeated integration by parts in the u variable 
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(gaining 2 ek 2 k l 2 from the non-stationary phase, but losing at most 2 ek 2 k / 2 from 
differentiating the <fi-(i/2-e)k function) then shows that this portion of the integral 
is 0(2~ ck ) for any C, and similarly for derivatives. Thus we may freely insert or 
remove a cutoff 4>-n/2-e)k( LJ ' ~ w ) m the integrand for W. In other words, we may 
restrict ourselves (smoothly) to the region 

uj = J + 0(2 £k 2- k / 2 ) (21) 

without any penalty. 

Now we argue similarly in the C variable, beginning with the identity 

Vq^x^AuXiV) =x-y. 
Thus if we insert the cutoff 1 — <j)_(i/2-e)k{ x ~~ v) m the integrand of W, we have 
\^c^x,k,z\ ^ 2( 1 / 2 ~ £ ) fc . Repeated integration by parts in the £ variable using (17) 
thus shows that this portion of the integral is 0(2~ ck ) for any constant C, and the 
same is clearly true for any £' derivatives of W. Thus we may restrict ourselves 
(smoothly) to the region^ 

y = x + 0(2 ek 2- k/2 ) (22) 

without any penalty. 

Now let us restrict ourselves smoothly to the regions @ and (||). Then we 
have 

V„¥ x ,fc, a (u;, C, y) = ? ■ V a V f $(y, u) - C 

= £' • V x V^$(x, J) - C + 0(2 £fe 2 fe / 2 ) 
= V x $(a;,C') -C + 0(2 £fc 2 fc / 2 ) 

Thus if we insert the cutoff 1 — <j>i 1/2+25) k(V x <&(:r, £,') — C) m the integrand of W, we 
have I ^y^ x ,k,z(i^, Ci y)\ ^ 2( 1 / 2+2e ' fe . Repeated integration parts in the y variable 
as before (noting that the gain 2~( 1 / 2+2e ) fc from the non-stationary phase exceeds 
the loss 2( 1 / 2_£ )' c coming from differentiating <j)-{i/2-e)k{ x ~ y)) thus shows that 
we can restrict ourselves smoothly to the region 

C = V x $(x, £') + 0(2( 1 /2+2e)fe ) (23) 

once we have first restricted to the regions ( ^l|) and (|2^). 
Because of the above restrictions, we have 

s(x, = s(x, V a $(x, O) + 0(2 Cefe 2- fc / 2 ). 

If we estimate the contribution of the error 0{2 Cek 2~ k l 2 ) by taking absolute values 
everywhere and taking full advantage of the restrictions ( ^l| ) , (^2|) , J23| ) , we obtain 
a bound of 0{2 Cek 2~ k l 2 2~^ n ~^ k l 2 ) as desired, with each derivative in £' gaining 
an additional 2 C£fc 2~' c . Thus we may replace s(x, Q by s(a;, V x Q(x, £'). 

Now we remove the </>(i/2+2e)fc(V a; < I>(:r, £') — () cutoff by reversing the above 
argument^. We can then evaluate the £ integration using the Fourier inversion 

6 An inspection of this argument shows that one can in fact localize further, to the region 
y = x + 0(2 ek 2~ k ). However if one does so, the uncertainty principle then prohibits one from 
localizing £ to any scale finer than 0(2~ Ek 2 k ), and one ends up with worse estimates at the end. 

7 Strictly speaking, this creates some formal difficulties because £ then ranges over all of R n 
and the oscillatory integral is not absolutely convergent. However one can fix this by applying a 
suitably large cutoff <£jc(C) an d eventually letting K —* +00, observing that with this cutoff the 
Fourier inversion formula holds in the weak limit. We ignore this technicality. 
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formula 

J e 2M*-y)< = S (y - x ) 
to simplify the above expression to 

s(x, V x $(x, £')) / e 2 ^'-^*^^)-^*^'-'))^^, w ) e -n{w)*m 

(1 - (p- £k (J(x, u)))\ J[x, w)| 1/2 0_ (1/2 _ e)fe (w' - cu) duj. 
From (|l5|), we observe the Taylor series expansion 



£' • (V{$(aj, w) - V { $(sb, w')) = ^A'V^ $(1, uj'){lu - w', « - w') + 0(2 Ce 2- fe/2 ) 

on the support of 0-(i/2- £ )/c( w ' — w )- The error term 0(2 Ce 2~ k l 2 ) will give us a 
contribution of 0{2 Cek 2~ {jl ~ X)k l 2 2~ k l 2 ) to VF (as can be easily checked by replacing 
every term by its absolute value) , with each derivative in £ picking up an additional 
factor of 0(2 Cek 2~ k ). Thus we may ignore this error. Similarly, we can use the 
Taylor expansion 

V {x,uj)e-^^/\l - <j>_ ek {J{x,u>)))\J{x,u)\ 1 / 2 
= pfowOe-'^V^tl - <t>- Ek {J{x, u/)))| J{x, oj')\ 1/2 + 0{2 Ce 2- k ' 2 ) 

(noting that fi must be constant, since J cannot change sign from uj to uj 1 ) and 
discard the error term 0(2 Ce 2 - *' 2 ) as before. Thus we may replace W with 

S (x,V^(x,C0)^,^0e^ VW4 (l-^(J(a ; ^')))l^^')r /2 

One can replace the cutoff function <^_(i/2- e )fc(w' — uj) by a more standard bump 
function 4>o(uj' ~ u) if desired, again by the principle of non-stationary phase. Stan- 
dard stationary phase asymptotics (see M) then give 

e w<A'^»(a,u;')(«-c-',a»-«-')0 o ( £( ;'_ w ) duj = e"' 1 '"''"/ 4 1 J(x, fc/) |~ 1/2 (A') -("-Wa+O^" 

for any C, and similarly for all derivatives in Comparing this against Wo we 



thus obtain (20) as desired. ■ 
The proof of Theorem [O] is now complete. 

7. Remarks 

It seems that there is some room for improvement in Theorem First of all, 
one can relax the angular regularity of the symbol a{x 1 £), and measure a in a more 
exotic symbol class. For instance, a natural class (which has appeared elsewhere, 
e.g. ||) seems to be those symbols obeying the estimates 

l^afojafooi < c QA7 ir^- |/3|+¥ - 

Also, when the phase function becomes degenerate (i.e., when J becomes small) 
the estimates above improve, roughly by a factor of \ J\ X / 2 , when £ is large. Of 
course, in the most extreme case, when $ vanishes, then the Fourier integral op- 
erator collapses to a pseudo-differential operator, and one has weak-type (1, 1) for 
operators of order 0, and not just — (n — l)/2. In intermediate cases when V 2 $ 
consistently has rank strictly between and n — 1, there are intermediate results 
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(see 0j, ||). Perhaps one can unify these results by introducing symbol classes 
adapted to the phase function <fr(x, £); for instance one might study symbols a(x, £) 
which obeyed bounds such as 

\a(x,0\ <(l + 2^ k \J(x,0\r 1/2 
together with some corresponding bounds on higher derivatives. It is not clear to 
the author exactly what the best symbol classes to use are, but the techniques here 
are likely to miss the optimal class by at least an e in the exponents. 
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